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In this paper we quantise scalar perturbations in a Randall-Sundrum-type model of inflation where 
the inflaton field is confined to a single brane embedded in five-dimensional anti-de Sitter space- 
time. In the high energy regime, small-scale inflaton fluctuations are strongly coupled to metric 
perturbations in the bulk and gravitational back-reaction has a dramatic effect on the behaviour of 
inflaton perturbations on sub-horizon scales. This is in contrast to the standard four-dimensional 
result where gravitational back-reaction can be neglected on small scales. Nevertheless, this does 
not give rise to significant particle production, and the correction to the power spectrum of the 
curvature perturbations on super-horizon scales is shown to be suppressed by a slow-roll parameter. 
We calculate the complete first order slow-roll corrections to the spectrum of primordial curvature 
perturbations. 

PACS numbers: 04.50.+h, ll.10.Kk, ll.10.St 

I. INTRODUCTION 

Recent developments have suggested that our four-dimensional Universe could lie on a brane embedded in a higher- 
dimensional space-time This new perspective may dramatically change our picture of the early Universe. One 
of the most studied scenarios is the model proposed by Randall and Sundrum where a single brane is embedded in 
five-dimensional anti-de Sitter (adS) space-time Q. In this model, the Friedmann equation is modified at high energy 
densities Q and the bulk gravity is strongly coupled to the dynamics on the brane. Because of that inflation in 
the early Universe, and the perturbations generated during inflation, could be significantly modified compared with 
conventional four-dimensional models. 

The simplest way to realize inflation in the brane world is to consider inflation driven by the potential energy of a 
scalar field, the inflaton, confined to the brane 0]. The amplitude of the resulting scalar curvature perturbations is 
given by 

(K 2 C ) 1/2 = (j) (WY /2 , (1-1) 

where H is the Hubble parameter, <f> is the time derivative of the inflaton </>, and 8(f) is the inflaton fluctuation on a 
spatially flat hypcrsurface. The quantum expectation value of the inflaton fluctuations on super-horizon scales in the 
de Sitter space-time is 

» j >=(!) 2 . M 

It has been shown that the constancy of the curvature perturbation TZ C on comoving hypcrsurfaces outside the horizon 
is independent of the gravitational theory Q and thus it is also valid in this brane- world model Q. The curvature 
perturbation, Eq. (|l.ip , can thus be directly related to the observables like Cosmic Microwave Background (CMB) 
temperature anisotropy. Although the formula for the curvature perturbation is exactly the same as in the four- 
dimensional case, the relation between the Hubble parameter H and the scalar field potential V is modified at high 
energy densities due to the modification of the Friedmann equation Q ■ This changes the prediction of the spectrum 
of scalar perturbations for a given potential. Based on Eq. (jl.lj) , interesting progress has been made on observational 
predictions of the brane- world inflation models @, [1] ■ 

A crucial assumption in deriving Eq. (|1.2p . and hence calculating the amplitude of the curvature perturbation (jl.ip . 
is that back-reaction due to metric perturbations in the bulk can be neglected. In the extreme slow-roll limit, where the 
coupling between inflaton fluctuations and metric perturbations vanishes, this assumption must be valid, but it is not 
obvious whether this assumption is still valid once slow-roll effects are taken into account. Ref. @ initiated the study 
of bulk metric perturbations generated by inflaton fluctuations on a brane. It was shown that sub-horizon inflaton 
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fluctuations on a brane excite an infinite ladder of Kaluza-Klein modes of the bulk metric perturbations already at 
the first order in slow-roll parameters. Subsequently Ref. (To| found that, due to this infinite ladder of modes, a 
naive slow-roll expansion breaks down in the high energy regime once one takes into account the back-reaction of the 
bulk metric perturbations. This was confirmed in Ref. pd| by direct numerical simulations that solved the coupled 
equations for inflaton perturbations and bulk metric perturbations without using any slow-roll approximation. 

Although these results indicate that the back-reaction of the bulk metric perturbations give rise to corrections of 
order one to the behaviour of inflaton fluctuations on sub-horizon scales even in slow-roll inflation, this does not 
necessarily imply that the amplitude of the inflaton perturbations receives corrections of order one on large scales. 
In order to address the amplitude of the inflaton perturbations, and hence the curvature perturbations (jl.ip on large 
scales, we need to quantise properly the coupled brane inflaton fluctuations and bulk metric perturbations. The 
purpose of this paper is to perform quantisation of the coupled fluctuations and determine the corrections to the 
inflaton fluctuations amplitude, Eq. (|1.2[) . In this way we study whether the bulk metric perturbations give rise to 
a significant modification of the amplitude of the primordial fluctuations, Eq. (jl.ip . generated during high-energy 
inflation on a brane. 

This paper is organised as follows. In section II, we begin by describing the model and present the coupled equations 
for the inflaton fluctuations on a brane and the metric perturbations in the bulk. Slow-roll parameters are defined 
and our approximation for dealing with the branc-bulk coupling is explained. Based on this approximation, late-time 
solutions are derived. In section III, we study the behaviour of the bound states, which determine the final amplitude 
of the curvature perturbations. We perform numerical simulations to evolve the late-time bound states backwards in 
time and to find their early-time behaviour. We also develop a semi-analytic method to derive solutions for the bound 
states at early times. The bound states are shown to consist of an infinite ladder of the Kaluza-Klein modes. Using 
the semi-analytic analysis, asymptotic solutions for perturbations on small scales are derived. Section IV is devoted 
to the quantisation of the bound states. We present a method to determine the quantum amplitude of fluctuations 
with the knowledge of asymptotic solutions for the bound states on small scales. Then the corrections to the formula 
(jl.ip at the first order in the slow-roll parameters are presented and compared with the standard slow-roll corrections 
in four-dimensional gravity. We end up with our conclusions in Section V. 



II. DESCRIPTION OF THE MODEL 



A. Basic equations 

We take as a starting point the analyses of Refs. 0, [T^, EH- The model under consideration is the single-brane 
Randall-Sundrum model @ with a cosmological constant in the bulk and an inflaton field <p confined to the brane, 
whose potential V(<p) satsifies the slow-roll conditions. The action for the model is 



(2.1) 



where K is the trace of the extrinsic curvature of the boundary 3A4 and A is the tension. The adS curvature scale is 
defined by A = — 6/i 2 and the tension is tuned as A = Qji/ k 2 . 
The background line clement has the form 

ds 2 = -N 2 (y,t)dt 2 + A 2 {y 1 t)8 lJ dx l dx j +dy 2 . (2.2) 

Here the warp function A(y, t) and the lapse function N (y, t) are Q 

A(y, t) = a(t) jcosh^y) - (l + ^ sinh(//y)| , (2.3) 
N(y, t) = cosh( M y) - jl - g(2p + 3p) J smh( M y), (2.4) 

p=^4? + V(<j>), p^^-V^). (2.5) 

In these coordinates, the brane is located at y = yi> = 0. The scale factor is determined by the Fricdmann equation 
and the equation of motion for the scalar field, 

2 

# 2 = yp(l + ^), 4> + ZH4> + V'{4>)=Q. (2.6) 



where 
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Allowing for the linearised metric perturbations gives the metric 



-N(y, tf{\ + 2A)dt 2 + A(y, t) 2 {l + 2Tl)6 ij dx i dx j + (1 + 2A yy )dy 2 + N{y, t)A y dydt. 



(2.7) 



It was shown in Refs fl2l [HI ] that the perturbed five-dimensional Einstein equations are solved in the adS background 
if the metric perturbations arc derived from a master variable O: 
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The remaining perturbed five-dimensional Einstein equation then gives the wave equation in the bulk, 
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(2.9) 
(2.10) 

(2-11) 
(2.12) 



for perturbations with comoving wavenumbcr k on the brane. 

On the brane, there is also the perturbation of the inflaton 8<t>. It is useful to express the scalar field perturbations 
Sep on the brane in terms of the gauge invariant, Mukhanov-Sasaki variable [l4j 



Q 



where subscript b denotes the quantity evaluated at the brane. The equation of motion is [111 ] 

2 
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+ V'\4>) } Q = J(O), 



(2.13) 



(2.14) 



where 



H = (A'/A) h and k\ = k 2 ^. Here 
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(2.15) 

(2.16) 

(2.17) 
(2.18) 



come form the projected Weyl tensor [H, [TB|, while the curvature perturbation, Eq. (|2.11[) . evaluated at the brane is 
given in terms of O as 



n h = — ( met' - 3hci - 3ri 2 ci + ^rCi 

6a 1 



(2.19) 



The junction condition, which describes how the brane field provides the boundary condition for the bulk metric 
perturbation is [ill ] 
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(2.20) 



The master equation (|2.12p . the Mukhanov-Sasaki equation (|2.14p with a source (|2.15p and the junction condition 
(I2.20[) describe the coupled system of the inflaton and metric perturbations. 
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B. Slow-roll inflation 



The slow-roll parameters are denned as 



e 



H 




(2.21) 



In terms of the scalar potential, these parameters are given by [3] 




H/ /z 2> 1), the modification of the Fricdmann equation eases the condition for slow-roll inflation for a given potential. 

The amplitude of the scalar perturbations is determined by the curvature perturbations on a comoving hypersurfacc, 
which is conserved on large scales: 



At the zeroth order in the slow-roll parameters, the quantum expectation value of squared Q is (Q 2 ) = (H/2ir) 2 . 
Slow-roll dynamics of the inflaton field gives corrections to this result. One contribution comes from the evolution of 
Q on super- horizon scales. As is seen from Eq. f|2 . 14[) . deviations from the de Sitter space-time give an effective mass 
to Q and modify the behaviour of Q outside the horizon. In four-dimensional cosmology this correction is known as 
the Stewart-Lyth correction. This Stewart-Lyth correction in the present brane-world inflation model was calculated 
in Ref. 0] assuming that the coupling to the bulk metric perturbations could be neglected. 

In four-dimensional cosmology, the slow-roll corrections do not affect the dynamics of Q on sub-horizon scales, as 
the effective mass induced by the dynamics of the inflaton can be neglected in comparison with the gradient term. It 
is then possible to quantise Q as an effectively masslcss field in the de Sitter space-time. However, in the brane-world 
case there is no reason to believe that the coupling to the bulk metric perturbations can be neglected on sub-horizon 
scales. In fact in Refs. [13, [H| it is shown that the coupling to the bulk metric perturbations can lead to 0(1) 
effect on the behaviour of Q on small scales at high energies (H/fi S> 1) even though the coupling is suppressed by 
the slow-roll parameter. However, it has been an open question how this coupling affects the quantum amplitude 
of Q on larger scales. The aim of this paper is to investigate the effect which the mixing between the bulk metric 
perturbations and the brane inflaton fluctuations has on the quantum amplitude of Q after crossing out the horizon. 
We arc only interested in corrections of the first order in the slow-roll parameters. Thus, we can study separately the 
effect of the brane-bulk coupling and other slow-roll corrections. Almost to the end of this paper we neglect all the 
slow- roll corrections in the equations of motion except for the coupling between Q and f2, as we study the effect of 
the brane-bulk coupling. 



As explained above, we assume zeroth-order slow-roll for the background everywhere except for the terms inducing 
the mixing between the bulk and brane perturbations. This means that, in practice, we consider the background as 
the de Sitter brane configuration [l7] ]. 




(2.24) 



C. de Sitter approximation 




(2.25) 



where 



N(z) 
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(2.26) 



(ismh(Hz) 



and the brane is located at 




(2.27) 
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Here, H is the Hubble parameter (constant in the de Sitter Universe) and r is the conformal time coordinate. Hereafter 
we redefine the time variable 



T —> —T . 



so that evolution forward in r is evolution backwards in time. 

It is convenient to define a new variable Xi related to the master variable as follows, 

n = k- 2 a 2 N 3/2 X - 
The equation of motion for \ derived from Eq. (|2.12p is 



(2.28) 



(2.29) 



where 



X,TT + k 2 x-^X+j^(~x" + U{z) X ) =0, 



U(z) = ^(9~— 1 



4 V sinh 2 (Hz) 

Note that Eq. (|2.30|) is separable. The other equations are simplified by defining a new brane variable 



(2.30) 



(2.31) 



u = aQ+ qh*\ 



It obeys the following equation of motion, 
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derived from Eq. ()2.14|) . and the junction condition 

X \ 2 N 3 



n\4>H [tu_ t — u] 



(2.32) 



(2.33) 



(2.34) 



derived from Eq. (|2.20p . Here we have neglected all terms suppressed by the slow-roll parameters except for the terms 
responsible for the coupling between u and Xj i- e -> the terms that induce the mixing between the brane and bulk 
modes. The latter terms at the first non-trivial order in slow-roll are 0((3), where 
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To the leading order in the slow-roll parameters, /? 2 can be written as 

-1/2 



(2.35) 



(2.36) 



This is the parameter that controls the strength of the coupling between the inflaton perturbation on the brane and 
the gravitational perturbations in the bulk. 

The second order action for the variables f2 and Q has been derived in Ref. In terms of u and x, the action is 
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From this action, the coupled equations (|2.30p . (|2 . 33[) and (|2.34|) can be derived. It is useful to define the time- 
independent Wronskian of two solutions, 



W(T,f) - «*V - uu* T - g|-(«*x| b - «X*| b ) +7^j™dz( 



X X,t XX, r J 



(2.38) 
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where / = [u, x) is a state comprising a brane and a bulk mode. The Wronskian will play an important role in 
quantisation. 

Without mixing, u and x would behave as canonical scalar fields in four dimensions and five dimensions, respectively. 
The existence of the mixing terms in the field equations and in the Wronskian makes the analysis complicated. We 
simplify the equations further by defining yet another new variable to describe the perturbations of the inflaton field 
on the brane, 



c = u t H — u y , , 

^ ' r 3#t b ' 



(2.39) 



which is related to the density perturbation on a comoving hypersurface [13[ . The brane equation of motion (|2.33p is 
now written as 



and the junction condition (|2.34|) becomes 



irr + fc 2 C = - 
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(2.40) 



(2.41) 



Upon substituting (|2.39p into ()2.38|) . the Wronskian becomes 



W(f*J) 



dz (x*X, 



XX, T ) 



(2.42) 



5 Jz 



where / = (£,%)• 

In this and next sections we are going to study the coupled system of equations (|2.30[) . (|2.33|) and (|2.34[) . or, 
equivalently, (j2"^D|) and (gUJ). 



D. Late time solutions 



At late times, it is possible to find analytic solutions to the coupled equations. Solutions divide into bound states 
and continuum modes. There are two independent bound state solutions @. A growing bound state solution is 

+ X, = C, T ^h(^)lo g ( ^g);; ), (2.43) 

where A g is the overall amplitude, while B and C g can be determined through A g . In terms of £, the growing mode 
solution is given by 

i 9 = k 2 A g , (2.44) 
which is clearly the growing mode solution of Eq. (|2.40[) at late times fcr — > 0. A decaying bound state is given by 



u d = A d r 2 , Xd = C d t 2 v /sinh(i?z) 

where Cd is linearly related to We will re-derive these solutions later on in a slightly different language. For the 
time being we do not need explicit relations between the coefficients entering (|2.43[) or (|2.45[) . 
On the other hand, the continuum modes are 

ux=g X T a \ Xx = fx(z)T a \ (2.46) 

where 

a A = i±zA, (2.47) 

and in the high energy regime 



cosh(_ff;z) / cosh(Hz) — 1 
1 + 2 ° g I cosh(Hz) + 1 



(2.45) 



h{z)=N x cos(\z + p x ). 



(2.48) 
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The coupled equations determine g\ and (3\ in terms of N\. Again, we do not need the explicit relations here and 
in what follows. It suffices to say that using these solutions, one can check that the Wronskian between the bound 
states and the continuum modes vanishes [23| as r — > 0, and hence at all r. In other words, they are orthogonal to 
each other. Furthermore, unlike the growing bound state solution, the continuum modes decay away as r — > 0, so 
they do not have any effect at super-horizon scales. Thus, we focus on the bound states in the rest of this paper. 

We are interested in the final amplitude of perturbations, which is determined by the the growing bound state. On 
large scales, the growing mode gives 

tt2 

TZ C = ^A g . (2.49) 

The amplitude A g has to be determined by quantisation of the inflaton and metric fields on sub-horizon scales. We 
thus have to evolve the growing bound state backwards in time. For reasons that will become clear later, the early-time 
behaviour of the decaying bound state is also of importance. 



III. BOUND-STATE SOLUTIONS 



A. Numerical analysis 



We now investigate what becomes of the bound states at early times. By using the numerical method developed in 
Ref. [HI, the coupled equations are evolved backwards in time using the late time bound state solutions, Eqs. (|2.43|) 
or (|2.45j) . as initial conditions. In this section, we rescale to dimensionless variables 



k~ 



H~ x z. 



(3.1) 



to simplify the problem further. Recall that the evolution backwards in time is the evolution forward in r because we 
redefined the time variable r — > — r. The boundary equations <|2.40[) and (|2.41| become 
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In terms of rescaled £ and \, the Wronskian is 
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w(tj), w(rj) = w c (r,o + w x (x*,x), 



where 



(3.3) 



w ( (C,0 = SXt - eCr. w x(x*,x) = W 2 / dz (x*x,r - xx: T ) • 



(3.4) 



Fig. 1 shows the behaviour of x(z,t) for the mode that grows as r — > 0, i.e., has late-time asymptot 
clear that the bound state persists at large r. As can be seen from Fig. 2, the amplitude of x 011 the 
for large r like x{ z bi T ) x t 1 ^ 3 . As for the profile in the bulk, x( z i T ) is more and more concentrated 
at earlier times (larger r). This means that the solution for %(z,t) is not a separable function of r 
other hand, the amplitude of £ approaches a constant for large r as is seen from Fig. 3. The right 
shows that the phase difference between the growing mode and the decaying mode is roughly equal 
shows the behaviour of the Wronskian between the growing mode and the decaying mode. For large r 
is conserved independently for and W x . In these calculations, we have set A g = Aj, = 1. 
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B. Kaluza— Klein tower 

The numerical results indicate that even at early times, there exists the bound states. In Refs. 0, [T^], the behaviour 
of the bound states has been analyzed by performing the pcrturbativc expansion in 1 . It was shown that, in terms 
of the Kaluza-Klein decomposition, the bound state is composed of an infinite ladder of modes with negative mass 
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FIG. 1: The growing mode solution for x{ z i T )- We take f3 — 0.1 and H/fi — 100. The brane is located at Z), — 5.3. 
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FIG. 2: The growing mode solution for x( z , T ) f° r the same parameters as in Fig. 1. Left: the profile of x( z , T ) m the bulk. 
The amplitude is normalised so that the amplitude on the brane is equal to 1. Right: behaviour of \ on the brane, x( 2 f>i r )- 



squared. Due to the excitation of the infinite ladder of the modes, the naive perturbation expansion in 1 breaks 
down at early times [Toj |. In this paper, we derive the solutions without performing the perturbative expansion in 1 . 

The bulk equation of motion for \, Eq. (|2 . 30[) is straightforward to solve. The solutions can be written in the 
form 



X = y dv B u f u {z)T u {T) , 



where v is related to the conventional Kaluza Klein mass as 



H 2 ' 
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(3.6) 



t _1 ' 2 T„(t) is in general a linear combination of the Bessel function J v (t) and Neumann function N v (t), and 
(sinh z)^ 1 / 2 f u {z) is a linear combination of the Legendre functions P„_i/2(ooshz) and Q v - 1/2 (cosh z). We are inter- 
ested in the bound states for which \ vanishes as r — > and/or z — > 00, so we consider the solutions composed of the 
following mode functions: 



T v {r)=T X ' 2 J v {r), f v (z) 



( sinhzX 1 ^ 2 (5 J ,_i/ 2 (coshz) 
Vsinhz b / Q l/ _ 1/2 (coshz b ) 



(3.7) 
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FIG. 3: Left: The growing mode solution for £(r) for the same parameters as in Fig. 1. Right: Comparison between the 
growing mode £ s (solid line) and decaying mode £<j (dotted line). 
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FIG. 4: Behaviour of the Wronskian between the growing and decaying bound states, W(f g , fd)- 



where Re v > 0. Our normalisation ensures f v {z b ) = 1- We can then write the brane field £ in terms of the mode 
functions T v (t) as 



£= / dvCMr). 



(3-8) 



The brane equation of motion (|2.40[) and the junction condition (|2.41|) then reduce to two coupled difference 
equations relating the coefficients B v and C„ entering Eqs. (|3.5[) and (|3.8[) . 
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where we use the shorthand 
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We now derive two independent bound state solutions to these coupled equations. 



(3.11) 
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To get a bound state we require that B v = whenever Re v < 0. It can be seen from Eqs. (|3.9p and (|3.10[) that 
there are two such solutions: one for which C 2rl _i/2 and B 2n+1 / 2 are non- vanishing only, and the other with non-zero 
values only of C2 n +i/2 and B 2n+3 / 2 , where n = 0,1,2, . ... Importantly, we choose C_i/2 real, then all C2„_i/2 and 
B 2n +i/2 obtained from Eqs. (|3.9p and (|3 . 1 0[) are real; the same is true for the second solution. 

At late times, the modes with higher n decay faster, thus the late time solutions are described by the lowest mode: 



Xg = B 1/2 f 1/2 {z)T 1/2 {T), £ g = C_ 1/2 T_ 1/2 (t), 
Xd = B 3/2 f 3/2 (z)T 3/2 (T), £ d = C 1/2 T 1/2 (t) , 



where C_x/2 and C\/ 2 arc arbitrary constants and according to 
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We see that the solution that starts from C_x/2 corresponds to the growing mode (|2.43|) with A 

the solution that starts from C\/ 2 corresponds to the decaying mode (|2.45p with A c i = Cx^v^/ 71 "- Hence, solving the 
recurrence equations (|3.9[) and (|3 . 1 0[) starting from C_\i 2 or C\i 2 is equivalent to evolving the late time bound state 
solutions backwards to obtain their early time behaviour. We have solved the recurrence relations numerically, the 
results are shown in Figs. 5 and 6. Using these solutions, we have checked that the solutions x( T , z ), £( T ) obtained 
by summing up the ladder of modes are in excellent agreement with the results of numerical simulations described in 
section IIII Al 
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FIG. 5: Solutions for C v with f3 2 — 0.1. The growing mode solution with 6*2,1-1/2 is shown by + and the decaying mode 
solution with C2 n +i/2 is shown by x. We took C1/2 = CL1/2 = 1. 




FIG. 6: Solutions for 6*2,1-1/2 an d B 2n +i/2 f° r large n. 
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C. High energy limit: analytical treatment 



From now on we will work in the high energy limit H/ fi 1. In this limit one has 

as H/fj, — > oo . 



and the expression (|3.7[) for the mode function simplifies to 



,-v(z-Zb) 



(3.14) 



(3.15) 



In order to derive analytically the asymptotic behaviour of the bound state solutions at early times, r — > oo, it is 
useful to have an asymptotic expansion for B v and C v for large n or cquivalently for large v. From the numerical 
solutions we find that the sequences -B 2n +i/ 2 and C 2rl _i/2 alternate in sign (and the same for the other solution) so 
we define new sequences 



i?2n+l/2 — (~l)™ +1 ^2n+l/2 i Cin-X/I — {~ l)"Can— 1/2 ? 

for the growing bound state and 

-^271+3/2 = (~ l)" +1 -32n+3/2 , C 2n+ i/ 2 = ( — l)™C2n+l/2 i 

for the decaying bound state. For large v, we can write the difference equations as 

dv 
d 



C v 



a: -1/2 
2x 



H r 



B„ 



-1/4 
2x 



X — V — l 
-■ X — V-\-\ 



C x 



dp 



iyc v ) , 



(3.16) 
(3.17) 

(3.18) 
(3.19) 



provided that B v and C„ are smooth functions of v whose derivatives vanish as v — ► oo. From Eqs. (|3.18[) and (|3.19D 
it follows that the coefficients B v satisfy the differential equation 



dB, 



2(5 2 B V = 0, 



dv V dv 

to the leading order in v^ 1 , which has solutions 

B v oc Z (2/?V2^) , 
where Zq is a linear combination of the Bessel function of order 0, 

Zo, i (2pV2i)=Aj, i J o (20V2Z)+A If ,iNo(20V2Z), 



(3.20) 



(3.21) 



(3.22) 



where i = g for the growing mode solution and i = d for the decaying mode solution, and Aji and A^^ are yet 
undetermined real constants, two for the growing mode and two for the decaying mode. The reality of Ajj and Ajq.i 
follows from the reality of B v . From (|3. 9|) we find, again to the leading order in v , 



(3.23) 



where Z\ is the combination of J\ and N\ with the same coefficients Aj^ and A^^ as in Eq. (|3.22p . In the large- n 
asymptotics, the coefficients for the growing bound-state solution are thus 



iW/ 2 = (-ir+^w^) , c 2n _ 1/2 = (-ir +i /3 



In the same way, the decaying bound state has 



£ 2n+3 /2 = (-l)" +1 ^o(4/3^) , C 2n+1/2 = (-l)» +1 /3 Zl(4 ^ , 



(3.24) 



(3.25) 
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We have checked that the agreement between these asymptotic solutions and the numerical solutions of Eqs. (|3.9|) 
and (|3.10[) is excellent after we fit the two coefficients Aji and A^,i for each mode. 

To the zeroth order in /3 2 , the constants Ajj and Ana can be expressed analytically through C_i/2 (for the growing 
mode) and C\i 2 (for the decaying mode). Since for the growing mode C„ is 0(/3 2 ) for v > —1/2, we can neglect C v 
with v > —1/2 in Eq. (|3. 9[) for the this mode. Likewise, we can neglect C„ with v > 1/2 for the decaying mode. 
Then, as was first shown in [§], Eq. Ij3.9| can be solved as 

iWi/2 = (-1)" Un + I) — — C_ 1/2 , B 2n+3 / 2 = (-1)" Un + |) — —C 1/2 . (3.26) 

\ "V ^2n+l/2 V ^2ri+3/2 

Making use of the expression (|3.14[) valid in the high energy regime, we get for n 3> 1 

B 2n+ i/2 = (-l)' l+1 2C_ 1/2 , B 2n+3/2 = (-l)' l+1 2C 1/2 . (3.27) 

Comparing with the non-perturbative solutions (|3 . 25[1 and (|3.24p . we find that the perturbative solutions (|3.2T|) are 
valid for /3y/n <g; 1. Taking the small argument limit of Eqs. (|3.25j) and (|3.24|) . we obtain 

A. lg = 2C_ 1/2 , A J)d = 2C 1/2 , A N , a = A Nd = 0, (3.28) 

at the zeroth order in /3 2 . The latter relations effectively provide matching between the late-time behaviour (|3. 121) 
and the early-time asymptotics which we are about to find. This matching, however, is valid to the zeroth order in 
1 only. 

It should be emphasized that the perturbative approach breaks down for (3y/n > 1. This is because, even if each 
C v with v ^ ±1/2 is 0{(3 2 ), the contribution to B v from C„ in Eq. ([33)1 is 0(n x C v ) ~ 0(nf3 2 ). As we will see 
later, large-n modes become dominant for large r. Then even for small /3 2 , the naive perturbation theory in j3 breaks 
down at large r. This is in accord with Ref. [1 01 ] . For finite /3 2 , Aj^ and A^.j are found from numerical solutions of 
Eqs. (|3.9p and (|3.10[1 and they are shown in Fig. 7. 

2.3 I 1 1 1 1 1 0.35 I 1 1 1 1 1 

2.25 - x X X X x x 0.3 - 



2.2 - X Aj,d 0-25 - 

2.15 x - 0.2 - 

2.1 - x + , - 0.15 - 

+ + 
+ + 

2.05 - + 0.1 - 

+ 

A, + * 

2S^ ■ rt - J >S -.. 0.05 - 

1.95 ' 1 1 1 1 ' 0*— 

0.02 0.04 0.06 0.08 0.1 



IA N ,d I 
IAw I 2 



I Aks I 

IAi, g l : 



(3 2 P 2 

FIG. 7: Numerical solutions for the coefficients Aj y i and |j4.jv,i| 2 /|^j,i| 2 , denned in Eq. (|3.22p . for different values of (3 2 . We 
took Ci/2 = C_i/ 2 = 1. 

Using the asymptotic expressions for the coefficients given in Eqs. (|3 . 25[) and (|3.24[) . we can approximate the 
solutions at z — Zf, <C 1 by 

Xg (r,z) ~ ^(-l) n+1 ^o(4/3V^)^2„ + i/2(r)e- 2 ™( z -^), 



t,z) ~ ^(-l)" +1 Zo(4/3V^)V^J 2 „ + 3/2(r) e - 2n(z -^ ) . (3.29) 



For large r, the sums over the modes can be evaluated explicitly. First we use the asymptotic expression for the 
Bessel functions, 



JJt) ~ \l 2 —— cosfwtanfl-i/0-7r/4), (3.30) 

>' -KV tan 6 x ' 
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where cos(9 = v/t, sec Rcf. [19fl . The sums are dominated by the modes for which n/r <C 1 but n 2 /t 3> 1. In this 
regime the asymptotic expression simplifies to 

/ 2 

Jn+i/2{r) = \ — cos (t + n 2 /(2r) - nn/2 - tt/2) . (3.31) 

V 7TT 

The standard large-argument asymptotic form [l9j for the Bessel function Zq in Eq. (|3.25p enables one to approximate 
the sums in Eqs. (pT2"§|) and (|3~jg)) as 



Xg oc (Aj^g + iAiy,g) J dnn' 1 ^ cxp (— 2n(z — z b ) — 4i/3\/n + it + 2in 2 r~ 1 ) + c.c. , 

Xd oc (A/ id + Mjv.d) J dnn^ 1 /* cxp(^-2n(z- z b )-4i/3y/n + iT + 2iii 2 T~ 1 -^i^J +c.c., (3.32) 

where we keep only the terms that give dominant contributions to the integrals (see below). The phase difference of 
7r/2 between the growing mode and decaying mode solutions comes from the fact that the growing mode and decaying 
mode solutions contain J2 n +i/2 an d ^2n+3/2j respectively. At large r (early times) these integrals are saturated by 
modes with large n. This means, in particular, that we cannot truncate the ladder of modes if we wish to describe 
the early time behaviour (Io| . The leading-order contribution comes from the saddle point at 

" saddlo = ^P V3 t 2/3 - |(* - z b )r . (3.33) 
The integrals in Eq. (|3.32|) can thus be approximated as 

Xg oc (A. lg + l A N , g y/z cxp (ir ^/? 4/ V/ 3 - ^^^(z - z b ) + l -r{z z h f} + c.c. , 

Xd oc {A. Ld + l A N . d )r 1 / 3 cxp (ir ^P A/Z r^ - 2 1 / 3 /3 2 /V 2 / 3 (z - z b ) + l -r{z - z b f ^ij + c.c. . (3.34) 

Note that at large r, these wave functions are localized in a narrow region in the bulk, z — z b ~ (/3t)~ 3 / 2 < 1, as we 
anticipated when writing Eq. (|3.29p . The terms omitted in Eq. (|3.32[) do not have saddle points near the semi-axis 
of positive n and thus give the contributions which are exponentially suppressed as compared to Eq. (|3.34p . In the 
same way we can find the asymptotic behaviour of £, 

£, g oc (Aj i3 + iA N>g ) cxp (it - J^/3 4/ V /3 ^ + c.c , 

Zd oc {A J>d + iA N . d ) cxp Ur - ^/? 4/3 r 1/3 - + c.c . (3.35) 

These asymptotic expressions nicely explain the behaviour of solutions at early times (t 3> 1) obtained in our 
numerical simulations. We see that the amplitude of x{ z bi T ) increases as x( z b, T ) °c t 1 / 3 while the amplitude of £ 
tends to a constant at large r. At the same time, x(z,t) is increasingly dominated by large-n modes and becomes 
more and more closely bound to the brane at large r. This behaviour ensures that the Wronskian for Xi given in 
Eq. (|3.4p . is conserved at late times. We also notice that the phase difference between the growing mode and decaying 
mode solutions is approximately equal to tt/2, provided that Aj^ 3> Ajy^. As we saw above, the latter inequality is 
indeed valid at small 1 , so the pattern shown in the right panel of Fig. 3 also has its analytical explanation. 

To conclude this section, we notice that at early times, the original coupled system Eqs. (|2.30p . (|2.40p and (|2.4ip may 
be treated within the WKB approximation. This is done in Appendix A. The WKB analysis is in prefect agreement 
with the results Eqs. I|3~34"|t and (|3~3"S]) . 



IV. THE QUANTUM AMPLITUDE 



In the preceding section, we have found the early time behaviour of the bound state solutions, Eqs. (|3.34[) and 
(|3.35p . For finite (3 2 they are parameterized by the real coefficients Aj^ and An,% where i = g,d for the growing and 
decaying bound state solution, respectively. These coefficients can be found numerically, see Fig. 7. Importantly, the 
expressions (|3.34[) and (|3.35p show that despite the complicated behaviour of the bound states at early times, the 
division of solutions into positive- and negative- frequency parts is well defined as r — > oo. This is also clear from the 
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WKB analysis, sec Appendix A. Thus, as in the familiar four-dimensional case, we can straightforwardly quantise the 
perturbations and define the vacuum state on small scales (early times) to determine the final amplitude of the inflaton 
fluctuations. We follow the method developed in Refs [lo, Hl[ to relate the amplitude of the growing mode solution 
to the creation and annihilation operators. In principle, this method involves the construction of properly normalised 
positive- and negative-frequency solutions at early times and the evaluation of the inner product (Wronskian) between 
these solutions and the bound states. However, there is a short-cut: all we need are the values of the real coefficients 
Aj t i and An,i- We first explain this short-cut by making use of the familiar four-dimensional example, and then 
proceed to the brane-world model. In this section, we restore k, but we are still working with redefined time r — > — r. 



A. Four-dimensional case 



Let us begin with the standard four-dimensional case. At the first order in the slow-roll parameters, the equation 
for u is 

u TT + k 2 u- -^(2 + 6e„ - 3r)„)u = 0. (4.1) 
In the same way as in Eq. (|2.39[) . we define 

it 



£ = u T + (l + 2e H - Vh )-. (4.2) 



T 



The equation for £ is 



The Wronskian is defined as 



trr + k 2 t- 2 Vh £ = 0, (4.3) 



;) = P«*6t-£C-)- ( 4 - 4 ) 



Neglecting k 2 in Eq. (|4.3[) . one obtains the growing mode and the decaying mode solutions at late times, r — > 0, 

L -> {kT) 1 ' 2 ~ l/ 2 v ^^- -, £j -> (kT) x ' 2+v 2- u ^ -. (4.5) 

where v = 1/2 + 2e H — r\ H and we included numerical factors for later convenience. We are interested in quantum 
field which behaves at late times as 

l^H a , (4-6) 

where Z is a time-independent quantum operator. This field is quantised on small scales where we can expand £ in 
terms of negative- and positive-frequency modes, 

1 = acp^ +oV (+) , (4.7) 

where a and arc annihilation and creation operators, respectively, which define the vacuum 

d|0)=0. (4.8) 

We should note that due to the redefinition of time r — > — r, the definitions of the negative- and positive-frequency 
functions are <p(~> oc e lkr and oc e~ lkT . The mode functions should be normalised as 

W(ip { -\<p {+) ) = -i . (4.9) 

The latter condition ensures the canonical commutational relation between £ and its conjugate momentum. 

Our aim is to express Z in terms of d and aJ. We exploit the constancy of the Wronskian for this purpose. At late 
times, Z can be related to the Wronskians, 
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Here W(£ g ,£d) can be calculated using the late time solutions (|4.5[) for £ g and i d 

1 /2 



W(e 9 ,^) = ^(-JsinK). (4.11) 

On the other hand, W(i,id) can be calculated at early times. For this purpose, we expand the growing mode and 
the decaying mode solutions in and tf\~> at early times, 

t g = c g ^ + c> (+) , £d = + c^ (+) . (4.12) 

Then, W(i,i d ) is calculated as 

W(iU) = -i(ac d -tfc d ). (4.13) 

As usual, the operator Z corresponds to the Gaussian random field, fully characterised by the vacuum expectation 
value of its square. The latter is evaluated as 

=fc 2 (!) 2 M 2 (sin^)- 2 . (4-14) 

Thus the problem reduces to finding c d , which is the expansion coefficient for the decaying mode solution. 

Now, the trick is that we do not need to know the precise expressions for to determine \c d \. Indeed, by 

solving for the growing mode and decaying mode backwards, it is possible, at least in principle (and in practice in the 
four-dimensional example), to derive the early-time solutions for £ d . In the four-dimensional case they are 

C 9 = N g cos(kr + A g ), £ d = N d cos(kT + A d ). (4.15) 

This immediately impics that the positive- and negative-frequency modes have the form 

p(±) = \ v \e* i5 e* ikr , (4.16) 

where the normalisation factor \ip\ may be found from (|4.9p . but will not be needed in what follows, and S is some 
phase (which is irrelevant in the four-dimensional case). Comparing Eqs. (|4.12[) and (|4. 1 5[) . we get 

^j = ^, Aigc g -Argc d = 5 g -6 d = A g -A d . (4.17) 
\ c g\ iv s 

The additional information on c d , c g and \5 d — S g \ can be obtained from the Wronskian W(i g , £, d )- In terms of c g and 
Cd, the Wronskian is calculated as 

W(C g ,U) = 2\c g \\c d \\sm(5 g -5 d )\. (4.18) 
Now we can use Eq. (|4TTj) together with Eqs (|4~T7)) and ((418)) to obtain \c d \ 2 in terms of N d /N g and |A a - A d |, 

Cd 2 = __£__ _ __ _ _ sini/7r. 4.19) 

1 d| 2N g |sin(A g - A d )| k \ir ) v 1 

Thus, the quantity \cd\ 2 entering (|4.14p is expressed in terms of the ratio Nd/N g of the amplitudes and difference 
(A g — Ad) of the phases of the decaying and growing modes evolved back in time. 
In the four-dimensional case, the growing mode and the decaying mode solutions are 

^ = (fcr) 1 / 2 J_„(fcr), i d = ikr) 1 ' 2 J v (kr), (4.20) 

At late times, r — > 0, they have the form (|4.5p . while at early times, fcr 3> 1, the solutions become 



^V: cos ( fcr + -v-K - 7 tt ] , i d -> \j - cos ( kr - -vtt - -it ) . (4.21) 

We thus find 



Nd 



1, A g - A d = isir. (4.22) 
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Then \cd\ 2 is obtained as 

1 (1 



M 2 = ^j, (4-23) 
and the expectation value for Z squared is given by 

(Z f Z) = ^(sini/7r)" 2 . (4.24) 

In this way, from the solutions for the growing and decaying modes, we can determine the quantum amplitude of the 
fluctuations. We note in parenthesis that the slow-roll corrections do not affect the value of \cd \ in the four-dimensional 
case. This is because the slow-roll effects do not play any role at early times in Eq. (|4.3p . 
The quantity that is related to observables is the comoving curvature perturbation, 

TZ C = rU. (4.25) 

a<j) 

Using the result for the quantum expectation value for Z< Z, Eq. (|4.24[) . and the late-time asymptotics (|4.5p . we get 
from (j4~2|) at fcr « 1 

(u 2 )V 2 _ 2 2, H -r lH F(3/2 + 26, Zjkl J= {kT) -^ H + n H . (4 . 26) 

Expanding the Gamma function in Eq. (|4.26|) . we find 

(u 2 ) 1/2 - {l + (2e H - r lH )(2 - 7 - ln2)} 7 L=(fcr)- 1 - 2 ^+^ . (4.27) 



'2k 

Here we have used the formula for the psi function 

V(3/2)=^|^=2- 7 -21n2, (4.28) 
where 7 is the Euler constant. The power spectrum of the curvature perturbation is then given by 

,*.(*) _£<*.V£(»)V). (4.29) 

Using the fact that conformal time is given, up to the first order in slow-roll parameters, by 

r = ^-(l + e H ), (4.30) 
all 

one finds the power spectrum of the curvature perturbation, 



H\ 2 ( iJ N 



l + 2S(e H ,r, H ) , S(e H , VH ) = (2e H -r, H )(2-ln2-j)-e H , (4.31) 



2tt / 

where S is known as the Stewart-Lyth correction [22j . 

B. Brane- world case 

We now apply the method presented in the preceding subsection to the brane-world case. The late time behaviour 
of £ is the same as in the four-dimensional case: at late times we have 

i^Z£ g . (4.32) 

We expand the quantum operators £ and x a t early times as 

I = a<f { - ] +aV + \ 

X = a<£ (-) +a t (+) , (4.33) 
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where and <j>^ are the positive and negative frequency solutions to the coupled equations, obeying the normal- 
isation condition 



(4.34) 



where = (tp^, 4>^), and the Wronskian has precisely the form (|2.42p . This normalisation ensures the canonical 
quantisation conditions for the fields, as is shown in Appendix B. 

The growing mode solution and the decaying mode solution are also expanded as 



6 



(4.35) 



where i = g,d. At late times, the contribution from x to the Wronskian can be neglected and the dominant contribution 
comes from the lowest mode in the ladder of the modes, Eq. (|3.12[) . 



£„ = (fcr)J_ 1/2 (fcr), £ d = (fcr)J 1/2 (fcT), 
where we set C\/2 = C-1/2 = 1- 111 this way we obtain the Wronskian, 

»u, /<) = -■ 

Then, as in the four-dimensional case, the expectation value of the quantum operator Z^ Z is calculated as 



(4.36) 



(4.37) 



(4.38) 



The virtue of the approach outlined above is that it works even in cases when the early-time solutions have complicated 
form. In the brane-world case, the positive- and negative-frequency solutions can be read off from Eqs. (|3.34|) and 
(|3.35p . and they indeed are rather complicated. Nevertheless, the above argument goes through, and we find, in exact 
analogy with the result (|4. 19)) , 



l_Nd 1 1 

2N g |sin(A s -A rf )| k 



(4.39) 



where 



A 2 



A 2 



A rl = - 



arctan ■ 



.4 



N,g 



A 



J,a 



arctan ■ 



.4 



N,d 



A 



j.d 



(4.40) 



As a cross check, let us show that the standard four-dimensional result is recovered at the zeroth order in f3 2 . At this 
order, the solutions for Aj yi and Anj are given by Eq. (|3.28|) . Then Eq. (|4.40p yields 



\N g \ 2 



= 1, An - Ad 



7T 

2' 



(4.41) 



and we get 



(4.42) 



which agrees with the four-dimensional result, Eq. (|4.23|) . 

For finite 1 we use numerical solutions for Aj^ and A^,i shown in Fig. 7 to evaluate Nd/N g and (A„ — Ad). The 
results are presented in Fig. 8. Using these numerical results, we obtain the expectation value of interest, 



1 + K{j3 2 ) 



(4.43) 



where K(f3 2 ) is shown in Fig. 9. As is seen from Fig. 9, K(f} 2 ) ~ 0(/3 2 ) ~ 0(e H ). 
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FIG. 8: Numerical solutions for \N d \/\N g \ and |A S - A d - ir/2\. 
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FIG. 9: Numerical solutions for K(f3 2 



C. Amplitude of curvature perturbations 



Finally, the amplitude of the comoving curvature perturbations is calculated using Eqs. (|2.49p and (|2.44p . At late 
times, £ behaves as 



i = z 



Then, comparing this with Eq. (|2.44[) , we find the amplitude of A 



Using Eq. (|2.44p , the spectrum of the curvature perturbations is evaluated as 



-, 9 , ^ 9 

\ ( n 

2^ 



1 + K(P 2 ) 



(4.44) 



(4.45) 



(4.46) 



We should note that we did not take into account the standard slow-roll corrections at the super-horizon scales 
in this calculation. The correction K{(3 2 ) is solely due to the coupling between the inflaton and the bulk metric 
perturbations. Although the mode functions are significantly affected by this coupling at early times, the final 
correction to the quantum amplitude at late times is suppressed by the slow-roll parameter even in the high energy 
regime. 
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As the correction K{0 2 ) is of the same order as the usual slow-roll corrections, we should also take into account the 
corrections to the amplitude at the superhorizon scales. Since K(f3 2 ) is already the first-order correction, the slow-roll 
corrections to K(0 2 ) are of higher order. Then we just need to add the Stewart-Lyth correction. Indeed, neglecting 
the coupling to the bulk metric perturbations that is already first order in the slow-roll parameter, the equation for 
u is given by 



k 2 u- -1(2 



6e. 



3r] H )u = 0. 



(4.47) 



We find that this is exactly the same as in the four-dimensional case. Thus, the Stewart-Lyth correction is precisely 
the same as well; this was first shown in Ref. Q. Then at the first order in slow-roll, we get the final result 



l + K((3 2 ) + 2S(e H ,rj H ) 



(4.48) 



where S = (2e H — r) H ) (2 — In 2 — 7) — e H is the Stewart-Lyth correction. Fig. 10 shows the complete first order 
corrections in the brane-world model at high energies, H 3> /1, compared with the Stewart-Lyth corrections. 
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FIG. 10: First order corrections 2S(e H ,t) h ) +K(f3 2 ) in the high energy regime ff >/i, compared with 2S(e H , r\ H ) for r\ H = 0.1 
(left) and r] H = 0.01 (right). 



V. CONCLUSIONS 



In this paper, we performed the quantisation of the inflaton perturbations during slow-roll inflation driven by an 
inflaton confined to a branc in the Randall-Sundrum brane world-model. The slow-roll dynamics leads to a finite 
coupling between the inflaton fluctuations on the brane and metric perturbations in the bulk. We first studied the 
effect of this coupling at the first order in the slow-roll expansion by keeping the coupling term while neglecting all 
other slow-roll corrections. This enabled us to deal with five-dimensional perturbations around the de Sitter branc 
for which the bulk wave equation is separable and solutions can be described as sums of bulk modes. At late times 
there exists bound states with single bulk modes, and the growing bound state determines the final amplitude of the 
comoving curvature perturbation on large scales. The continuum modes were shown to be orthogonal to the bound 
state in the sense that the Wronskian between them vanishes; furthermore, the continuum modes decay away at late 
times. Thus we could neglect the continuum modes in our analysis. 

We evolved this late time bound state backwards numerically to find its early time behaviour. It was shown that 
the behaviour of the inflaton fluctuations is significantly affected by its coupling to an infinite tower of bulk modes, 
even at the first order in the slow-roll approximation. Asymptotically, the amplitude of £, that describes the inflaton 
perturbations on the brane, becomes constant while the amplitude of \, that describes the bulk metric perturbations, 
increases on the brane at early times. It was also found that x is more and more concentrated near the brane at 
early times, becoming dominated by modes with large (and negative) Kaluza-Klein mass-squared. We found that it 
is possible to reproduce these numerical solutions by the summation of the infinite ladders of modes. The growing 
mode and decaying mode solutions have different ladders of modes, which leads to a phase shift between them. The 
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infinite sum of the modes can be evaluated using a saddle point method and analytic formula for the asymptotic 
solutions was derived. 

We then performed the quantisation of the bound state. We developed a scheme where the quantum amplitude of 
the perturbations can be derived by finding the early time behaviour of the growing mode and the decaying mode. 
In this method, the constancy of the Wronskian is exploited to connect the quantum operator at late times to the 
creation and annihilation operators defined at early times, without the explicit analysis of the early-time positive- 
and negative-frequency functions and the calculation of the Bogoliubov coefficients. This method was applied to the 
standard four-dimensional case to rc-derive the Stcwart-Lyth correction, which arises from corrections to the field 
evolution after horizon exit. We then applied this method to our brane-world model. In this case, there is a correction 
that arises from the modified solutions for the perturbations on small scales. Thus, this correction is due to the 
coupling of the inflaton to the bulk metric perturbations. Despite the fact that the behaviour of perturbations on 
small scales is qualitatively different due to this coupling, it was found that particle production is suppressed by the 
smallness of the coupling. The correction to the final amplitude was thus shown to be of the order of the slow-roll 
parameter, e H , even in the high energy regime. 

Because the effect of the mixing between the brane inflaton fluctuations and the bulk metric perturbations is of 
the first order in the slow-roll parameter, the complete first order corrections consist of the correction K(f3 2 ) coming 
from the branc-bulk mixing and the Stewart-Lyth correction S(e H , f] H ), which was shown to have the same form as in 
the four-dimensional case. Fig. 10 is our main result showing the complete first order corrections in the brane-world 
model. 

In four-dimensional general relativity the corrections to the evolution of inflaton fluctuations that arise at the first 
order in the slow-roll parameters can be neglected for wavelengths much smaller than the Hubble scale. Thus, the 
vacuum fluctuations on sub-horizon scales are essentially the same as those of a massless field in the dc Sitter space- 



time 14J. The only change in the amplitude of quantum fluctuations comes from the modification of evolution on 



super- horizon scales, leading to the familiar Stewart-Lyth correction [221 ] . 

By contrast, in the brane world the coupling between the inflaton fluctuations and bulk metric perturbations that 
arises at the first order in the slow-roll parameters leads to a radically different picture of the bound sates at early 
times compared with the zeroth-order solution. Short wavelength inflaton fluctuations do not decouple in the high 
energy regime (H 3> /it) from the bulk metric perturbations and support an infinite tower of bulk modes, which 
become increasingly dominated by largc-n modes on small scales. Nonetheless, the mixing between positive- and 
negative-frequency modes is absent in the sub-horizon regime, while additional mixing near the horizon crossing is 
suppressed by slow-roll parameters, so the final correction to the amplitude of fluctuations on super- horizon scales is 
again first order in the slow-roll parameters. 
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APPENDIX A: WKB APPROXIMATION FOR EARLY-TIME SOLUTION 



Let us consider early times, i.e., large r. In this case we can find the bound state solutions to the coupled system of 
equations (|2.30[> . (|2.40|) and (|2.41[) by making use of the WKB approximation. Using the rescaling given in Eq. (|3.1[) 
we approximate the solutions by 

e = A(r)ex P (<5c(r)) J x = B(t, z) exp(iS x (r, z)) , (Al) 

where A and B are slowly varying functions but the exponential factors are varying rapidly. Consistency with the 
brane equation of motion (|2.40[) requires that 



s d T ) = S x( T , z b) = S b (r). 



(A2) 
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In the WKB regime and in the high-energy limit, Eqs. (|2.40[) and (|2.41[) reduce to 



dS, 
dr 



1 A 



2j3 A 



dz 



B=tA. 



From these equations we obtain the boundary condition for S x , 



dS-y 



2i(3 2 



Let us also write the equation of motion for \, Eq. (|2.30[) . at the boundary, 



OS, 
dr 



1 + - 



9S V 



The latter two equations give 



dS, 
dr 



= 1 - 



2/3 2 



dz 



dS* 



= z2 1 / 3 /3 2 / 3 r 2 / 3 



One can use this approach to solve for S x as a Taylor scries about the brane position. The next derivative is 



dz 2 



(A3) 
(A4) 

(A5) 
(A6) 
(A7) 
(A8) 



and so on. Note that this solution indeed satisfies the applicability conditions for the WKB approximation. 

We can obtain the amplitudes A and B in a similar way by using the next order of the WKB approximation, finding 
that 



A = const , B. = const • r 1 / 3 

b 



(A9) 



Note that the amplitude of £ is constant but the amplitude of x increases like r 1 / 3 . The expansion near the brane 
gives 

X (t, z) oc rV3 e^ 1/3 " 2/V/3 (^) e -* b . (A10) 
This expression implies, in particular, that the contribution from the small- z region to the normalisation integral 

dx(x*X,r-c.c.), (All) 

is constant in r. The asymptotic solution Eq. (|A10j) agrees with the solution obtained in Section III.B, Eq. (|3.34j) . 



APPENDIX B: QUANTISATION CONDITION FOR COUPLED SYSTEM 

We normalised the positive- and negative-frequency modes as 

i i r°° 

We justify this as follows. In terms of £ and Xi the quadratic action is given by 



1 



1 



= ^2 / - k 2 e) + ^ / dzdr ( x 2 T - k 2 X 2 - 



v 5 Jz b 



X 12 + (U(z) - 2H 2 )x 2 
(Hr) 2 



1 1 fN' 



dr 



(Hr) 



(Bl) 



(B2) 
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Therefore, the conjugate momenta are 



H = *x = -^iX.r, (B3) 



where the minus sign comes from the redefinition of time r — > — r. From the expansions of the fields in terms of a 
and d', 

£ = atp ( ^ + aV (+) , 

X = a^ (_) + a f (+) , (B4) 

we can express a and a* as 

W{ft-\f) = -ia\ W{f^\f)=ia. (B5) 
Then the commutation relation between a and a' is calculated as 

[a, at] = [w(f( + \f),w(f(-\f)) 

+ p | d2(^+^(-) - ^ T "V +) )[4>x] + 3^2/ <te(#V _) " 4"V (+) )[X,^]. (B6) 
Imposing the canonical commutational relations 

[C,7rd=«! [x( z )> 7r x( z ')] = i^z,z'), [Ci 7r x]= > [X,i"d= ' ( B7 ) 
and using the normalisation (jBlj) we get 

[a,a f ] = l. (B8) 

This justifies our normalisation. 

A cross check is the calculation of the Hamiltonian which is derived from the action (|B2j) in the standard way. By 
making use of the decomposition (|B4|) . the Hamiltonian for the bound state is calculated as 



H k = k (aa f + , (B9) 

provided we normalise the modes as in (|B1|) . Thus, with our normalisation we obtain the free harmonic oscillator 
picture as we should. 
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